﻿FLAC 3D实例
FLAC 3D, 实例, 练习
1.网格的建立
new
gene zon brick p0 0 0 0 p1 40 0 0 &
p2 0 80 0 p3 0 0 30 size 10 10 6
model mohr
property bulk 0.5e8 shear 0.2e8 friction 40
property cohesion 1e10 tension 2e10
set gravity 0 0 -20.00
initial density 500
fix x range x -0.1 0.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 79.9 80.1
fix z range z -0.1 0.1
plot add surface yellow
plot set rotation 35 0 30
plot add axes red
plot show


2.监控最大不平衡力
new
gene zon brick p0 0 0 0 p1 40 0 0 &
p2 0 80 0 p3 0 0 30 size 10 10 6
model mohr
property bulk 0.5e8 shear 0.2e8 friction 40
property cohesion 1e10 tension 2e10
set gravity 0 0 -20.00
initial density 500
fix x range x -0.1 0.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 79.9 80.1
fix z range z -0.1 0.1
history nstep 10
history unbal
set mechanical force 100
solve
plot history 1
print history
plot show


3.监控某一点的位移变化
new
gene zon brick p0 0 0 0 p1 40 0 0 &
p2 0 80 0 p3 0 0 30 size 10 10 6
model mohr
property bulk 0.5e8 shear 0.2e8 friction 40
property cohesion 1e10 tension 2e10
set gravity 0 0 -20.00
initial density 500
fix x range x -0.1 0.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 79.9 80.1
fix z range z -0.1 0.1
history nstep 10
history gp zdisplacement 4 8 3
set mechanical force 100
solve
plot history 1
print history
plot show







4.位移等值线图
new
gene zon brick p0 0 0 0 p1 40 0 0 &
p2 0 80 0 p3 0 0 30 size 10 10 6
model mohr
property bulk 0.5e8 shear 0.2e8 friction 40
property cohesion 1e10 tension 2e10
set gravity 0 0 -20.00
initial density 500
fix x range x -0.1 0.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 79.9 80.1
fix z range z -0.1 0.1
history nstep 10
history unbal
history gp zdisplacement 4 8 3
set mechanical force 100
solve
plot create trench
plot add contour disp
plot set rotation 35 0 30
plot add axes black
plot show



5.应力等值线图
new
gene zon brick p0 0 0 0 p1 40 0 0 &
p2 0 80 0 p3 0 0 30 size 10 10 6
model mohr
property bulk 0.5e8 shear 0.2e8 friction 40
property cohesion 1e10 tension 2e10
set gravity 0 0 -20.00
initial density 500
fix x range x -0.1 0.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 79.9 80.1
fix z range z -0.1 0.1
history nstep 10
history unbal
history gp zdisplacement 4 8 3
set mechanical force 100
solve
plot create trench
plot add bcontour szz
plot set rotation 35 0 30
plot add axes black
plot show







6.任意剖面上的等值线图
new
gene zon brick p0 0 0 0 p1 40 0 0 &
p2 0 80 0 p3 0 0 30 size 10 10 6
model mohr
property bulk 0.5e8 shear 0.2e8 friction 40
property cohesion 1e10 tension 2e10
set gravity 0 0 -20.00
initial density 500
fix x range x -0.1 0.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 79.9 80.1
fix z range z -0.1 0.1
history nstep 10
history unbal
history gp zdisplacement 4 8 3
set mechanical force 100
solve
plot create trench
plot create GravV
plot set plane dip 90 dd 90 origin 4 8 3
plot add boundary behind
plot add bcontour szz plane
plot set rotation 35 0 30
plot add axes black
plot show







7.开挖后1000次迭代的位移等值线图
new
gene zon brick p0 0 0 0 p1 40 0 0 &
p2 0 80 0 p3 0 0 30 size 10 10 6
model mohr
property bulk 0.5e8 shear 0.2e8 friction 40
property cohesion 1e10 tension 2e10
set gravity 0 0 -20.00
initial density 500
fix x range x -0.1 0.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 79.9 80.1
fix z range z -0.1 0.1
history nstep 10
history gp zdisplacement 4 8 3
set mechanical force 100
solve
plot create trench
plot create GravV
plot set plane dip 90 dd 0 origin 8 16 6
plot set rotation 35 0 30
plot add boundary behind
plot add bcontour szz plane
save Trench.sav
plot current Trench
property cohesion 1e3 tension 1.5e3
model null range x 20 24 y 0 80 z 20 30
set large
initial xdisplacement 0 ydisplacement 0 zdisplacement 0
step 1000
plot create DispCont
plot copy GravV DispCout settings
plot add contour disp plane behind shade on
plot add axes
plot show







8.开挖后的应力分布图
new
gene zon brick p0 0 0 0 p1 40 0 0 &
p2 0 80 0 p3 0 0 30 size 10 10 6
model mohr
property bulk 0.5e8 shear 0.2e8 friction 40
property cohesion 1e10 tension 2e10
set gravity 0 0 -20.00
initial density 500
fix x range x -0.1 0.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 79.9 80.1
fix z range z -0.1 0.1
history nstep 10
history unbal
history gp zdisplacement 4 8 3
set mechanical force 100
solve
plot create trench
property cohesion 1e3 tension 1.5e3
model null range x 20 24 y 0 80 z 20 30
set large
initial xdisplacement 0 ydisplacement 0 zdisplacement 0
step 1000
plot add contour disp 
plot set rotation 35 0 30
plot add axes
plot show 
